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Description 

This invention relates to measuring multiple optical 
properties of biological specimens, such as a population 
of cells, e.g., blood cells, at high rates of speed using 
computer controlled instruments that scan the biological 
specimens. 

Examples of such computer controlled instruments 
include flow cytometers, automated blood cell analyzers 
and blood cell differential classifiers. 

A flow cytometer is an instrument that hydrodynam- 
ically focuses a fluid suspension of cells into a single file 
stream that passes through an examination zone. A fo- 
cused light beam illuminates the cells in this zone and 
the instrument measures optical interactions of the light 
with the cells such as, for example, multiple wavelength 
absorption, scatter as a function of angle, and fluores- 
cence as a function of either wavelength or polarization. 
This type of instrument permits the study of living cells 
in addition to those which have been chemically treated, 
for example, by staining. Flow cytometry techniques 
certain constituents or structures.particularly those 
present on the cell surface, to be quantitatively charac- 
terized at cell rates of a thousand or more cells per sec- 
ond. 

Blood cell analyzers typically consist of a computer- 
ized microscope that automatically classifies various 
types of white blood cells and flags and counts all ab- 
normal cells in a specimen. Using such an automated 
instrument, the operator can also view the cells manu- 
ally, stop the analyzer temporarily for making visual mor- 
phological observations, or review abnormal cells in 
greater detail. During the automatic count, the X-Y co- 
ordinates of every encountered abnormal blood cell are 
stored. Therefore, the operator has the option of using 
a review mode in which abnormal cells are automatically 
and individually acquired and focused, and can view 
these cells on a TV. monitor or through a binocular mi- 
croscope. 

A blood cell differential classifier typically consists 
of a computer-controlled microscope having a stage that 
is driven by stepper motors. A light source, such as a 
xenon arc lamp, illuminates cells and the classifier uses 
various sensors such as a silicon photodiode array to 
measure the optical interaction between the cells and 
the light. All of the cells in a given area are illuminated 
at the same time. 

US Patent No: 4 647 531 describes a generalized 
cytometry instrument for measuring the kinetic proper- 
ties of live cells. A laser is employed to scan across a 
surface. The scan operates to find the position of a cell. 
Thereafter, the cell is illuminated, which may be by the 
same instrument, and an optical response or interaction 
is detected and recorded. Thereafter, the scan is re- 
sumed to find the position of a new ceil, and the process 
repeated. The spot size may be of a diameter commen- 
surate with that of the cell. 

We describe below a method and apparatus for 



generating optical data that accurately estimates multi- 
ple constituents and simultaneously characterizes a 
number of morphological properties of each of a popu- 
lation of cells. This method includes the steps of, and 
5 means for scanning the cell population with a beam to 
produce digital data samples, the different digital data 
samples representing multiple optical measurements at 
different locations within the cell population; storing the 
digital data, e.g., in a computer memory; locating a cell 
10 within the population, for example by comparing digital 
data derived from the stored digital data to a preselected 
threshold value; defining a neighborhood around the lo- 
cated cell; estimating a background level or individual 
background levels for all sample points in the neighbor- 
is hood based upon stored digital data corresponding to 
locations outside the neighborhood; and correcting 
each of the digital data samples corresponding to the 
neighborhood with the estimated neighborhood back- 
ground level to generate the optical data. The beam 
used in the invention is of electromagnetic radiation, e. 
g., laser light, X-rays, or infrared radiation. 

The sample correcting step may be performed by 
subtracting the corresponding background level from 
each of the stored digital data samples in a neighbor- 
hood. 

The term "digital data sample" describes digital in- 
formation that may be stored in a computer memory and 
that is sampled from an analog optical signal from the 
population of cells in a specimen when scanned with the 
beam. The optical data results from correction of the dig- 
ital data for background. Further processing of the opti- 
cal data generates corresponding optical and morpho- 
logical property values of the population of cells that may 
be used for precise measurements of cellular constitu- 
ents and as representations of cellular morphology. 

The term "constituent" as used herein refers to a 
specific part or component of a cell such as, for example, 
the DNA or RN A. Our methods are capable of determin- 
ing the amount of a specific constituent in a cell or cell 
population with an accuracy of within about a few per- 
cent of the actual amount. 

The term "neighborhood" means a specific region 
of sample points that is created by the operator to en- 
compass the digital data and that is set to statistically 
correspond to a single cell. That is, on average, only one 
cell will be located within a neighborhood. The size of 
such a neighborhood may be varied by the operator de- 
pending upon the size and concentration of the cells in 
the specimen. The neighborhood is used to allow an ac- 
curate estimate of cell constituents by determining the 
background within each neighborhood, i.e., individually 
for each cell. This allows for greater accuracy than a 
background level used for a number of different cells, 
because the background level varies across the exam- 
ining surface, such as a microscope slide. 

We also describe a method and apparatus for gen- 
erating optical data that characterizes a population of 
cells that includes the steps of, and means for, scanning 
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the cell population with a beam with a spot size compa- 
rable in size to the cells to be scanned to produce digital 
data samples, in which different samples represent dif- 
ferent locations within the cell population and the rate 
of sampling is such that the distance traveled per time 
as the spot is scanned between consecutive sampled 
locations approximates the spot size; and then process- 
ing the digital data to generate optical data that charac- 
terizes the cells in the population. The term "compara- 
ble' as used herein means within one order of magni- 
tude of the size of the cells to be scanned. 

In general, the number of sample points required to 
represent a cell is equal to the number of samples taken 
along a scan times the number of scans. When a small 
spot size of, e.g., one-half micron is used, as in the prior 
art, on the order of four hundred samples are required 
to accurately represent the constituents in a 10 micron 
by 10 micron cell area. We can now have a spot size 
which is on the order of 10 microns in diameter and is 
sampled at a rate such that the distance traveled per 
time as the spot is scanned between consecutive sam- 
pled locations is approximately the spot size. This gives 
the benefit of processing a sample about four hundred 
times faster than with the small spot size without losing 
accuracy. As the size is decreased and the sampling 
rate is increased this benefit correspondingly decreas- 
es. 

This method may also include the steps of multiply- 
ing the neighborhood value by a velocity normalization 
factor proportional to the velocity at which the beam is 
scanned past each location to generate optical data cal- 
ibrated for scanning velocity variations. 

We also describe a method and apparatus for gen- 
erating calibrated optical data that characterizes a pop- 
ulation of cells including the steps of, and means for, 
measuring the resulting intensity of a beam having a 
predetermined incident illumination intensity used to 
scan the cell population as a function of location of each 
sample within the scan; scanning the cell population 
with the beam to produce digital data samples, different 
samples representing different locations within the cell 
population; storing the digital data, e.g., in memory; nor- 
malizing the stored digital data using the beam intensity 
measurements of the location of each sample within the 
scan; and then processing the normalized digital data 
to generate calibrated optical data that characterizes the 
cells in the population. The resulting beam intensity may 
be measured, for example, with a sensor, such as a pho- 
tomurtiplier. 

We also describe a method and apparatus for gen- 
erating calibrated optical data that characterizes a pop- 
ulation of cells, that includes the steps of and means for, 
scanning the cell population on a surface with a beam 
to produce digital data samples, the beam having a giv- 
en illumination intensity, and different samples repre- 
senting different locations within the cell population; 
storing the digital data, e.g., in a memory; scanning a 
region of the surface without cells to generate intensity 



calibration data; generating an intensity normalization 
factor from the intensity calibration data and storing the 
intensity normalization factor in memory; locating a cell 
within the population, e.g., by comparing the digital data 
5 to a preselected threshold; defining a neighborhood 
around the located cell, the neighborhood containing 
optical data derived from the samples of digital data 
within the neighborhood; and correcting the optical data 
in a neighborhood, e.g., by multiplying the optical data 

10 by the intensity normalization factor, to generate the cal- 
ibrated optical data. 

The surface without cells may be a calibration slide 
including a uniform fluorescent dye surface, a dye-filled 
cuvette or a uniform scattering surface. 

'5 in another embodiment, the processing step may 
include multiplying the digital data by a velocity normal- 
ization factor proportional to the velocity at which the 
light beam is scanned by each location to generate op- 
tical data calibrated for scanning velocity variations. 

20 The method may also include the further steps of 
estimating background levels for the individual neigh- 
borhoods based upon stored digial data corresponding 
to each optical parameter at locations outside that 
neighborhood; and correcting each of the digital data 

25 samples corresponding to the neighborhood to gener- 
ate optical data, e.g., by subtracting the corresponding 
neighborhood background level from each of the sam- 
ples in the neighborhood. 

For each of the above methods, the neighborhood 

30 background levels can be determined by finding mini- 
mum digital data sample values in memory locations 
corresponding to actual locations adjacent the neigh- 
borhood for each optical parameter of each sample in 
the neighborhood, or by averaging a plurality of digital 

35 data sample values in memory outside the neighbor- 
hood. Typically, the term ■adjacent" means within five 
scans from the bounds of a neighborhood. 

In a further embodiment, time corresponding ap- 
proximately to the time the digital data was acquired is 

40 recorded and stored as synchronous digital time data, 
whereby each sample has a corresponding digital time 
point. As used herein, the terms "time corresponding ap- 
proximately to the time the digital data was acquired" 
mean the time a given strip of scans is completed. The 

45 cell population may also be located on a movable sur- 
face, the scanning may proceed in synchrony with a 
stepping motion of the surface, and the location corre- 
sponding to a sample may be derived from digital loca- 
tion data that may be generated, whereby each sample 

so has a corresponding digital time point and location. The 
movable surface may be a microscope slide, a cuvette 
into which the cells are inserted, or any other surface 
upon which or into which cells may be deposited. The 
cells on the movable surface of all of the above methods 

55 and apparatus can be visually observed. 

In another embodiment, an additional step of res- 
canning the population of cells, generating new digital 
data, location data and time data, and comparing the 
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new data with the corresponding data from a previous 
scan to determine changes in the cell population may 
be added. 

In a further embodiment of the above methods and 
apparatus, a cell can be located by detecting a digital 
data sample corresponding to the edge of the cell first 
contacted by the beam. In another embodiment, the 
neighborhood is initially centered around the first edge 
sample. The method may further include the steps of 
determining the maximum digital data sample within the 
neighborhood and then recentering the neighborhood 
with respect to the maximum sample before estimating 
a background level for the neighborhood. 

We also describe a method and apparatus for 
measuring multiple cellular properties including the 
steps of, and means for, scanning the cell population 
with a beam to produce sets of digital data samples, dif- 
ferent samples representing different locations within 
the cell population, e.g., storing the digital data in a 
memory; locating a cell within the population, e.g., by 
comparing the digital data to a set of preselected thresh- 
olds; defining a neighborhood around the located cell; 
summing all samples of optical data derived from the 
sets of digital data in a neighborhood to generate optical 
neighborhood values for the parameters for that neigh- 
borhood; comparing each parameter sample in a neigh- 
borhood with a set of preselected thresholds to generate 
a binary neighborhood pattern of samples either above 
or below the thresholds for each parameter; storing the 
optical neighborhood values and binary neighborhood 
patterns in memory; and generating optical properties 
corresponding to the optical neighborhood values and 
morphological properties corresponding to the binary 
neighborhood patterns for each cell. 

The optical properties may be generated by apply- 
ing math functions to one or more of the optical neigh- 
borhood values to generate optical property values of 
the cell corresponding to the neighborhood. The mor- 
phological properties may be generated by applying 
boolean functions to one or more of the binary neigh- 
borhood patterns to generate morphological property 
values of the cell corresponding to the neighborhood. 

The term "optical neighborhood value" describes 
the sum of all optical data in a given neighborhood for 
one parameter. This optical neighborhood value may 
then be used as is or processed further to generate op- 
tical property values proportional to specific cell constit- 
uents corresponding to cells of the population. The term 
"binary neighborhood pattern" describes a word pattern 
in memory corresponding to a single neighborhood, i.e., 
cell, and a specific parameter, that is generated by com- 
paring each sample point in a neighborhood to a thresh- 
old and registering a "1 " if the point is above a given 
threshold or a "0" if not. By testing each point, a memory 
word of ones and zeroes is generated which is the binary 
neighborhood pattern. This pattern can be further proc- 
essed to generate morphological property values of the 
cell. 



Each of these methods can include the further steps 
of estimating a set of background levels for the neigh- 
borhood based upon digital data corresponding to loca- 
tions outside that neighborhood; and subtracting the 
5 neighborhood background levels from each of the sam- 
ples of digital data corresponding to the neighborhood 
to generate optical data. 

We also describe an apparatus for measuring mul- 
tiple cellular properties that includes a beam source for 
10 generating a beam of electromagnetic radiation, e.g., 
light, X-rays or infrared radiation that may be controlled 
by the computer during the scanning process to illumi- 
nate the cells with a beam spot comparable in size to 
the cells to be scanned, an optical signal being gener- 
ated as a result of the spot illuminating the cells; a sur- 
face upon which the cells are located; an optical path 
for directing the spot from the beam source to the cells 
on the surface; a scanner that is interposed in the optical 
path between the light source and the cells for scanning 
the spot across the surface so that the spot passes over 
a plurality of cells on the surface; one or more sensors 
for measuring the optical signals; an analog-to-digital 
converter arranged to sample the optical signals at a 
specified rate and to produce digital data; whereby the 
sampling rate is set such that the distance traveled per 
time as the spot is scanned between consecutive sam- 
pled locations approximates the spot size; data proces- 
sors for correcting the digital data in the memory, gen- 
erating optical property values from corrected data 
points corresponding to multiple constituents of individ- 
ual cells, and processing the values to determine multi- 
ple cellular morphological properties; and a memory for 
storing the digital data produced by the converters. 

In further embodiments, the spot size used in the 
apparatus and in our methods is 3 to 50 microns or may 
be adjusted to be comparable in size to the size of the 
cells to be scanned. In addition, the wavelength distri- 
bution of the beam may be controlled by the data proc- 
essor. 

The scanner of the apparatus may be a resonant 
galvanometer scanner. Furthermore, the surface upon 
which the cells may be located may be a microscope 
slide, a cuvette into which the cells are inserted, or any 
other surface upon which or into which the cells may be 
deposited. 

In another embodiment of the above apparatus and 
methods, the celt population may be on a movable sur- 
face, the scanning may proceed in synchrony with a 
stepping motion of the surface and digital location data 
may be generated from which the location correspond- 
ing to a sample can be derived. 

The digital data produced by the cells and the beam 
in any of the above methods and apparatus may be de- 
rived from sets of analog optical signals. Furthermore, 
these optical signals may be fluorescence, forward an- 
gle light scatter, extinction, or wide angle light scatter. 

In another arrangement, we describe a method for 
generating calibrated optical data that characterizes a 
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population of cells in which a number of cells of that pop- 
ulation has a constant optical value for a given param- 
eter. This method corrects for beam illumination and 
light collection differences along each scan line. This 
method includes the steps of scanning the cell popula- 
tion with a beam along a scan line to produce digital data 
samples, different samples representing different loca- 
tions along the scan line within the cell population; stor- 
ing the digital data; locating a cell within this population; 
defining a neighborhood around the located cell, the 
neighborhood containing optical data derived from the 
stored digital data samples within the neighborhood; de- 
termining the location of the maximum digital data sam- 
ple within the neighborhood and then recentering the 
neighborhood with respect to the maximum sample; 
summing all optical data values in a recentered neigh- 
borhood to generate an optical neighborhood value for 
that neighborhood; determining the most frequently oc- 
curring optical neighborhood value to define a mode op- 
tical value; selecting a subpopulation of cells with optical 
neighborhood values within a predetermined range 
around the mode optical value; determining a set of av- 
erage optical neighborhood values of cells in the sub- 
population as a function of cell location along the scan 
line; computing an array of correction coefficients equal 
to the ratio of the mode optical value to the average op- 
tical value for each location along the scan line to gen- 
erate an illumination normalization factor; and correct- 
ing optical data for each cell at different locations along 
the scan line by multiplying the optical neighborhood 
value for that cell by the corresponding illumination nor- 
malization factor. The term "mode optical value" de- 
scribes the optical neighborhood value that occurs most 
frequently in a complete sample run. The term "illumi- 
nation normalization factor" describes an array of cor- 
rection coefficients equal to the ratio of the mode optical 
value to the average optical neighborhood value of all 
cells for each location along the scan line. 

In another arrangement, we describe a method for 
generating calibrated optical data that characterizes a 
population of cells. The method includes the steps of 
scanning a population of calibration particles comprising 
a known constant parameter with a beam along a scan 
line to produce digital calibration data samples, different 
calibration samples representing different locations 
along the scan line; storing the digital calibration data; 
locating a calibration particle within the population; de- 
fining a neighborhood around the located particle, the 
neighborhood containing optical calibration data de- 
rived from the stored digital calibration data samples 
within the neighborhood; determining the location of the 
maximum digital calibration data sample within the 
neighborhood and then recentering the neighborhood 
with respect to the maximum digital calibration data 
sample; summing all optical calibration data values in a 
recentered neighborhood to generate an optical neigh- 
borhood calibration value for that neighborhood; deter- 
mining the most frequently occurring optical calibration 



neighborhood value to define a mode optical calibration 
value; determining a set of average optical calibration 
neighborhood values as a function of particle location 
along the scan line; computing an array of correction co- 

s efficients equal to the ratio of the mode optical calibra- 
tion value to the average optical calibration value for 
each location along the scan line to generate an illumi- 
nation normalization factor; scanning the cell population 
with a beam along the scan line to produce digital data 

io samples, different samples representing different loca- 
tions along the scan line within the cell population; stor- 
ing the digital data; deriving optical data from the stored 
digital data; and correcting optical data for each cell at 
different locations along the scan line by multiplying the 

15 optical data for that cell by the corresponding illumina- 
tion normalization factor. 

Accordingly, there is provided, in accordance with 
a first specific aspect of the present invention, a method 
for generation cell constituent data that characterizes a 

20 population of cells, the method comprising: 

(a) scanning the cell population with a laser beam 
to produce digital data samples, different said sam- 
ples representing different locations within the cell 

25 population; 

(b) storing said digital data; 

(c) locating a cell within said population; 

(d) defining a neighborhood of digital data samples 
that includes the digital data samples correspond- 

30 jng to the located cell; 

(e) estimating a background level for each digital 
data sample of the neighborhood based upon 
stored digital data from corresponding scan loca- 
tions in prior or subsequent neighborhoods; 

35 (f) correcting each of said digital data samples cor- 
responding to the neighborhood with the estimated 
background level; and 

(g) summing said corrected digital data samples in 
the neighborhood to generate the cell constituent 
40 data. 

The invention provides, in a second specific aspect 
thereof, a method for characterizing a population of 
cells, the method comprising: 
45 (a) scanning the cell population on a surface along 
a scan line with a laser beam with a spot size compara- 
ble in size to the cells to be scanned, an optical signal 
being emitted as a result of said spot illuminating said 
cells; 

so 

a surface upon which the cells are located; 
an optical path for directing said spot from said 
beam source to the cells on said surface; 
a scanner interposed in said optical path between 
55 said beam source and the cells for scanning said 
spot across said surface along a scan line; 
a sensor for measuring said optical signal; 
an analog-to-digital converter arranged to sample 
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said optical signal at a specified rate and to produce 
digital data; whereby said sampling rate is set such 
that the distance travelled per time as the spot is 
scanned between consecutive sampled locations 
approximate said spot size; 

means for displacing the scan line of said laser 
with respect to said surface such that said laser 
spot scans continuously over a plurality of cells 
on said surface; 

a data processor for deriving cell constituent 
data from said digital data corresponding to in- 
dividual cells, and processing said cell constit- 
uent data to display cellular properties of the 
cell population; and 

a memory for storing said digital data produced 
by said converter. 

The invention is hereinafter more particularly de- 
scribed with reference to preferred embodiments there- 
of illustrated by the accompanying drawings. 

The drawings will first briefly be described. 

FIG. 1 is a schematic diagram of an instrument for 
measuring multiple optical properties of biological spec- 
imens at high rates of speed. 

FIG. 2 is a block diagram of the electromechanical 
circuit used in the instrument shown in Fig. 1. 

FIG. 3 is a flow chart of the general optical signal 
data acquisition loop. 

FIG. 4 is a schematic drawing of the scan pattern 
of the laser beam. 

FIGS. 5a and 5b are graphs illustrating the proper 
sampling rate according to the invention. 

FIG. 6 is a series of graphs representing analog op- 
tical data and the corresponding digital data. 

Fl G. 7 is a flow chart of the data processing function 
steps used to manipulate data stored in memory. 

FIG. 8 is a monitor screen display for an actual test 
run charting red fluorescence value (x-axis) against for- 
ward angle scatter value (y-axis). 

Fig. 9 is a flow chart of the data processing steps of 
a calibration method. 

Mechanical and Optical Systems 

Referring to Fig. 1, the instrument 10 includes a light 
source 12, a mirror scanner 14, such as a resonant gal- 
vanometer scanner, an epi-illumination microscope 16, 
a stepper motor controlled stage 18, light detectors 20, 
22, 24, 26, and various associated optical components 
which will be described below. Light source 1 2 produces 
a light beam 12a that reflects off of scanner 14 to pro- 
duce a scan beam 12b and finally illuminates a scan 
spot 1 2c of a fixed diameter or size on a specimen plane 
or surface 28. Specimen plane or surface 28 is posi- 
tioned on stage 18. Light source 12 is a laser such as, 
for example, a Helium Neon, Argon ion, Helium Cadmi- 
um, or solid state laser, depending on the application. 



More than one laser may be used for a given application. 
If this is the case, the beams can be combined using a 
dichroic mirror so that they are coaxial. For some appli- 
cations it may be desirable to control the intensity of the 
5 laser beams or shutter them under control of a compu- 
ter. Lasers with multiple wavelength outputs may also 
be used, in which case, it may be desirable to use a com- 
puter-controlled filter, prism or Bragg cell to select a spe- 
cific wavelength. 
10 After passing through a dichroic mirror 30, the laser 
beam is imaged by two lenses 32 and 34 onto an epi- 
illumination field stop 36 of microscope 16. Resonant 
scanner 14 is located between lenses 32 and 34 and 
scans the beam across the field stop when electrically 
75 driven. The focal lengths of lenses 32 and 34 and the 
deflection angle of scanner 14, which is proportional to 
the galvanometer drive voltage, control the size of the 
spot and the length of the scan at field stop 36 and thus, 
at specimen surface 28. The scanner is driven at 800 
Hz and the spot size at the specimen is 10 u. and the 
scan length at the specimen is 1 20 u,. These are nominal 
values and can be changed by the user by rotating mi- 
croscope nosepiece 38 bearing objectives 39 from the 
nominal 20X to other higher magnifications to reduce 
the spot size and scan length or lower magnifications to 
increase them. 

Epi-illumination is used to illuminate the specimen 
and to transmit fluorescent or scattered light to the view- 
ing eyepiece 40 or to a film or video camera 42. The light 
transferring assembly 44 may contain a dichroic or par- 
tially or fully silvered mirror as well as an optical filter in 
the viewing path. These assemblies are interchangea- 
ble and the microscope used in the described embodi- 
ment includes a movable rod to exchange these assem- 
blies. 

The specimen surface 28 may be a slide upon which 
a tissue or cytology specimen is mounted or a cuvette 
into which a cell specimen is injected. Furthermore, mi- 
croscope 16 may be an inverted microscope and spec- 
imen surface 28 may be a container or dish containing 
living cells. The optical instrument is designed to accept 
either a standard or inverted microscope. 

The cells of the specimen may be stained to alter 
their light scattering characteristics or cause them to flu- 
oresce. The scatter or fluorescence of each cell is relat- 
ed either to 1 ) the amounts of specific cellular constitu- 
ents such as DN A or proteins or 2) the amounts of one 
or more fluorescinated antibodies which bind only to 
specific cell antigens. 

The position of the specimen with respect to scan- 
ning beam 1 2b is controlled by X and Y stepping motors 
46, 47. These motors are driven in steps that are of a 
size ranging from a fraction of the spot size to the spot 
size. This step size is under the control of a computer 
program. The stage is also provided with sensors to pro- 
vide signals to the computer to indicate a "Home" or ref- 
erence position io for the stage and to limit its travel. By 
always moving the stage to Home at the beginning of 
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each run or periodically during a run (i.e. , a group of scan 
strips), it is possible to obtain the absolute position of 
given cells on the specimen surface so they can be man- 
ually reviewed or remeasured. For normal use, the 
stage can be moved manually by knobs attached to the 
motors. The stage may also be manually controlled by 
buttons, a joystick or a computer mouse with a specifi- 
cally programmed display. The microscope 16 is fo- 
cused manually by knob 62. 

Light scattered from the specimen in the forward di- 
rection is primarily related to cell size but may be altered 
by staining the cell. This scattered light passes through 
the microscope condenser lens 48, through a dichroic 
or partially silvered mirror 50, and is imaged onto two 
sensors 20 that are electrically connected but spaced to 
allow the incident scan beam image to pass between 
them. Therefore, only light scattered by the specimen in 
the forward direction away from the axial light path is 
collected by sensors 20. The microscope can be used 
normally to illuminate a specimen by the standard illu- 
minator 52 because the partially silvered or dichroic mir- 
ror 50, if suitably chosen, will not significantly reduce 
illumination. Mirror 50 may be placed into an assembly 
that mounts directly onto the microscope substage illu- 
minator and is small enough so that it does not interfere 
with focusing the condenser. The signals from sensors 
20 are amplified and become one input for the data ac- 
quisition circuit described in more detail below. Although 
not shown, an additional sensor may be placed behind 
or below sensors 20 to collect light from the incident 
scan beam in order to measure extinction. 

Light scattered at 90 degrees from the cells is useful 
for distinguishing certain types of cells such as blood 
granulocytes, because cells with granules scatter light 
at oblique angles. Obliquely scattered light will be 
trapped, with the specimen slide or cuvette acting as a 
light pipe, and can be detected by a sensor 22 placed 
at the slide or cuvette edge. This sensor is mounted on 
the stage in contact with the slide or cuvette and its am- 
plified signal is used as a second input for the data ac- 
quisition circuit. Again, this sensor does not interfere 
with normal microscope usage. 

Backscattered or 1 80 degree scattered light and flu- 
orescent light are collected by objective lens 38 at high 
numerical aperture and imaged back through the lenses 
35, 34, and 32 to dichroic mirror 54 which is designed 
to reflect the longer wavelength fluorescence if back- 
scatter is not to be measured or to reflect some of the 
laser source wavelength if backscatter is measured. Di- 
chroic mirror 30 transmits almost all of the laser wave- 
length. Dichroic mirror 54 splits the light into two parts 
so as to measure backscatter at the laser wavelength 
and longer wavelength fluorescence, or two different 
wavelengths of fluorescence. Mirror 54 can be partially 
silvered or reflect and transmit light based on polariza- 
tion if applications require measuring fluorescence de- 
polarization. Mirror 54 reflects part of the incident energy 
through filter 56 of the appropriate bandpass wave- 



length onto photomultiplier 24. Similarly, mirror 58 re- 
flects the remaining energy through filter 60 to photom- 
ultiplier 26. The signals from the two photomultipliers are 
amplified and become additional inputs for the data ac- 
5 quisition circuit. Apertures have not been shown in the 
figures and are used in various places to reduce stray 
light reflected from the surfaces of the lenses and mir- 
rors. 

10 Electromechanical System 

Referring to Fig. 2, the electromechanical system 
provides means of inputting signals from the light sen- 
sors 20, 22, 24 and 26, the scanner driver 70, and the 

'5 microscope stage 1 8, to the computer 81 and outputting 
signals from the computer to the X stepper motor 46 and 
Y stepper motor 47 which move the microscope stage 
18. This is accomplished by a commercially available 
circuit board 80 which accepts four analog voltages, dig- 

20 jtizes them at rates up to 100,000 Hz with analog-to- 
digital converter 80a and causes the resultant computer 
words to be stored in the computer memory under direct 
memory access (dma) control. Board 80 also accepts 
digital values from inputs 82-90, which provide limit in- 

25 formation from the stage, and provides digital output val- 
ues on lines 92-1 04, which control the stage. It also con- 
trols the values of two analog voltages used to control 
the supply voltages to photomultipliers (light detectors) 
24 and 26. A detailed description follows. 

30 In one embodiment, circuit board 80 is a Data 
Translation Model 2828 (Marlborough, New York) which 
plugs into one of the slots on an 80286 IBM compatible 
computer. The computer is represented in Fig. 2 by 
block 81. 

35 Furthermore, photomultipliers 24 and 26 are used 
to detect fluorescence or backscatter and photosensors 
20 and 22 are used to detect forward and oblique scatter 
emissions from the specimen and are each connected 
to respective amplifiers 76, 72, 78 and 74. These ampli- 

40 tiers have appropriate gains to provide signal levels of 
-1 0 to +1 0 volts to circuit board 80 at analog signal inputs 
106, 108, 110 and 112. The circuits 72 to 76 are D.C. 
operational or instrument amplifiers with D.C. "no sig- 
nal" levels near -10 volts. As will be described below, it 

45 is necessary to synchronize the position of the scanning 
mirror 1 4 with the data acquired by the sensors and con- 
verted to a stream of digital data. The digital data stream 
may be stored as two blocks of 64,000 words each. For 
synchronization purposes, the D.C. level of circuit 78 is 

so initially set to 0 volts, so that a negative synchronization 
pulse is easily detected. 

Synchronization is accomplished by using a square 
wave pulse synchronization signal generated by the 
scanner mirror driver 70 which controls the motion of 

55 scan mirror 1 4 through scanner 72. This signal is mod- 
ified to a pulse by circuit 114 and, in one embodiment, 
is added to the sensor signal from sensor 20. Although 
this reduces the digitization resolution of this sensor in- 
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put by two, a separate data channel for synchronization 
is not required. Of course, the synchronization signal 
may be used as a separate input signal. The signal at 
input 112 is the forward scatter signal of sensor 20 and 
the pulse signal is negatively added near one scan ex- 
treme (as shown in Fig. 6). As will be seen below, this 
negative pulse is detected by the program and used to 
properly synchronize the digital data stored in the com- 
puter memory. The sampling rate is set by the user 
through an initialization program which allows the user 
to define a Protocol for each test. The Protocol is a mon- 
itor screen that the operator uses to set the sampling 
rate and the various test parameters, area scanned, 
threshold settings, etc. The number of parameters dig- 
itized is also preset and the amplifier gain settings and 
input/output relationship, i.e., linear or logarithmic, may 
be used as additional parameters. 

The levels of the digital outputs 92, 94, 96, 98, 100, 
102 and 104, of circuit board 80 are under the control of 
the computer program. The digital inputs 82, 84, 86, 88, 
and 90 are read at specific times also determined by the 
control of the program. These outputs and inputs are 
used to control the movement of the microscope stage 
1 8 via X stepper motor 46 and Y stepper motor 47, which 
are each driven by translator circuits 116 and 117, re- 
spectively. The microscope stage is provided with limit 
switches which indicate when the stage has reached its 
limit of travel in the x and y directions. These switches 
generate signals on lines 18a-18d which are used as 
inputs 82 to 88, respectively, of board 80. Input 90 sens- 
es the pulses from scanner 70. Although not shown 
here, additional digital outputs may be used to control 
the wavelength of the light source by controlling specific 
light sources, shutters, or filter positions. 

The program controlled stage motion is designed to 
perform the following sequence that is depicted in the 
flow chart of Fig. 3. First, when the user initiates a test, 
both stepper motors are moved to a specific ■Home" po- 
sition. This is accomplished by calling a program sub- 
routine to pulse lines 94 and 1 00 on and off until inputs 
82 and 88 indicate that the stage has reached Home. 
An OR gate 118 passes the signal directly from output 
94 to input 1 1 6a. Under program control, outputs 96 and 
1 02 are set to produce the proper stage direction by pro- 
ducing signals received at inputs 116b and 117b. As 
soon as the stage reaches Home, the Y stepper is 
pulsed to move the stage 1 8 to the initial y position, then 
a subroutine is called to move the stage to the right in 
the x direction by changing the signal on output 96. In 
one embodiment, the pulse rate on output 94 is ramped- 
up in rate from about 100 up to 800 pulses per second 
(pps) for a fixed total number of pulses or distance, typ- 
ically 100 pulses or steps. This is the ramp-up number. 
This fixed ramp- up may also be adjusted by the program 
or by a commercially available ramp-up controller circuit 
(Metrobyte, Taunton, Massachusetts). 

The program produces a rate of 800 pps, which is 
the scanner frequency, at the end of the ramp-up, so the 



stepper is moving at full stepping speed at this time. The 
step size may be either 5 or 10 microns or some other 
size depending on the level at output 98 which controls 
the translator 116 through input 116c to produce full or 
s half steps of motor 46. Thereafter, output 92 is pro- 
grammed to provide an input to AND gate 120 to allow 
the stepper motor to be directly driven by the synchro- 
nization output of mirror scanner driver 70 stepping 
stage 18 to the right in the x direction at 800 steps per 
10 second, e.g., one step for each scan cycle. The stepping 
motion in conjunction with the scanning motion gener- 
ated by scanner driver 70 which is perpendicular to the 
stepping motion creates the scan pattern shown in Fig. 
4. In Fig. 4, the scan starts at the left, at the Home po- 
'5 sition, and forward and reverse scans of length "L" are 
produced until the end of one scan strip is reached. Such 
a strip typically encompasses 2500 forward and 2500 
reverse scans. An additional circuit may be provided be- 
tween driver 70 and translator 116 to allow other ratios 
20 of scanner rate to stepper rate under control of the user 
through the Protocol. 

An additional parameter of the Protocol is the X 
scan distance. This distance determines the length of 
one scan strip. This length can be used to calculate the 
25 size required for the data buffers by multiplying the 
number of parameters measured by the total number of 
data values digitized per scan strip. The board 80 digi- 
tizes inputs, stores them in a buffer, and returns a flag 
when the buffer is full. At this time, output 92 is set to 
30 transfer control of the X stepper from mirror driver 70 to 
program controlled output 94 via gates 1 20 and 1 1 8. The 
X stepper is then ramped down in velocity to a stop in a 
number of pulses or steps equal to the fixed ramp-up 
number, e.g., 100 steps. The digital data in the two buff- 
35 ers may be processed at the end of a complete scan 
strip, as will be described below, or it may be processed 
as it is being digitized. At the conclusion of this scan 
strip, stepper motor 47 moves the stage in the y direction 
so that a new scan strip can be run. Output 100 is used 
40 to send pulses to input 117a to step the motor; output 
102 determines whether movement is in the positive or 
negative y direction and output 104 determines the size 
of the steps (5 or 1 Ou.) and passes a signal to input 1 1 7c. 
Again, smaller or larger movements may be used under 
45 Protocol control. 

After stepping the Y motor 47 to move the specimen 
"up" or "down" in the y direction a distance equal to 60% 
of the scan length L, as shown in Fig. 4, the procedure 
described above, in which the X stepper motor 46 is 
50 ramped up in rate and moved the X distance, and 
ramped back down to a stop is repeated but the stage 
is moved back to the left in the x direction. The Y stepper 
then moves the stage in the y direction a number of 
steps and the complete cycle is repeated. The number 
55 m of Y steps are counted by the program and the test is 
terminated when the Y distance (a given M) is reached 
as determined by the user through the appropriate Pro- 
tocol parameter. 
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System Setup and Calibration 

Upon starting the program, the user is asked to se- 
lect a Protocol either from a set of named Protocols 
stored on disk memory, which can be modified at that 
time, or to generate a new Protocol. The Protocol is a 
set of parameters and demographic information about 
the tests to be performed on the sample. Demographic 
information describes the type of sample, the number 
or name of the sample, and any comments relating to 
the sample. The Protocol, along with the cell count for 
the test, is stored with each set of test results to properly 
identify the data and test conditions. The Protocol pa- 
rameters are: 1) the four start and end coordinates of 
the scan area and the X and Y step size, 2) the analog 
data acquisition rate, 3) the number of sample points 
processed along each scan line, 4) for each parameter, 
its name, its gain, a threshold level for cell findings a 
morphology threshold used to generate each cell's mor- 
phology word (binary neighborhood pattern), a flag to 
indicate which parameter will be used to find the cell's 
center, and where its values are displayed on the mon- 
itor during the test, 5) the size of the neighborhood 
(neighborhood is a concept to be described below), 6) 
the maximum cell count at which the test may be termi- 
nated, 7) the number of complete scans and 8) the prefix 
name of the file in which the test data is stored. 

Before beginning any sequence of tests, the user 
may be required to develop a set of calibration data by 
methods described in further detail below. The user can 
select calibration data stored from a previous calibration 
run or can generate new calibration data. If new calibra- 
tion data is desired, the user selects a section of the test 
specimen without cells or uses a specimen with no cells, 
places it on the microscope stage, and observes that 
the area to be scanned contains no cells or debris. The 
output of each parameter over the scan length is dis- 
played continuously as a function of time on the monitor 
screen during calibration and stage position can be 
manually controlled. To determine these background 
values, the user keys the beginning of a background cal- 
ibration scan and the values of 10 the parameters are 
digitized at one Y stepper position and one or more X 
stepper positions for a fixed number of scans on the or- 
der of 100. In one embodiment, this is accomplished in 
that the stage moves only a few steps in the x direction 
while continuously scanning up and down over gener- 
ally the same blank (no sample) area. The parameter 
values of each of the scan sample point positions are 
averaged over all scans. The average values and their 
coefficients of variation are displayed on the monitor 
screen to be checked by the operator to determine 
whether the machine is operating within the proper 
range of values. The coefficient of variation is the stand- 
ard deviation of a measurement divided by the mean of 
that measurement. These values are stored in memory 
as an array to be used in calculating cell parameter val- 
ues if required. 



The user can also initiate a second calibration run 
in which a test specimen provides an optical signal such 
as fluorescence or scatter. For example, a stained 
ground surface slide may be used. This is independently 
5 scanned as described above to derive a set of normal- 
ization factors for one or more parameters for each of 
the sample points. The calibration slide is scanned 
twice, the second time with the beam source blocked, 
so that the normalization factors are equal to the differ- 
to ence in the optical signal with and without the illuminat- 
ing beam. This eliminates D.C. and other signal noise 
in the system that is present even in the absence of il- 
lumination. If a given parameter is not normalized its fac- 
tor is one. These normalization factors are used to cor- 

15 rect the data values for variations in the optical signal 
received from the cells due to, e.g., the cell's location 
on the slide, the variation in the incident laser light, or 
other variations not related to a cell's constituents. 
These factors are then multiplied by the data points to 

20 remove the effect of these variations from the final 
measurements. 

The normalization factor may also be multiplied by 
a velocrty normalization factor to account for the differ- 
ent velocities along each scan because if most of the 

25 scan is utilized, more samples are digitized for cells near 
the ends of each scan length L because the beam slows, 
stops and turns in the opposite direction before starting 
the next scan length L The velocity calibration factor is 
the cosine of the distance of a sample point from the 

30 center of a scan length and is multiplied by the normal- 
ization factor for each sample point and stored as an 
array to be used later to normalize each sample point 
value in calculating cell parameter scan values. 

In another embodiment, an initial calibration of the 

35 device may be avoided by a method of self-calibration 
that operates during the sample run. This method avoids 
the need to scan a portion of a test specimen without 
cells or the use of any special calibration slide. The 
method generates an illumination normalization factor 

40 used to calibrate the optical data to correct for beam il- 
lumination and light collection differences along each 
scan line (length L) across the sample. When the cells 
of the sample are used to generate the normalization 
factor, the method is self-calibrating. In the alternative, 

45 however, the method can be used with calibration par- 
ticles in an initial run in which the calibration particles 
are scanned for a constant optical value of one of a 
number of known constant parameters, such as size or 
flurescence, to generate the illumination normalization 

50 factor. 

This method may be used when the cell population 
of the sample includes a number of cells that have a 
constant optical value for some parameter, such as, for 
example, the io DNA value of cells in the resting state 
55 or the size of the cells. To achieve an accuracy of about 
three percent in the calibrated data, approximately one 
thousand cells of the population having this constant pa- 
rameter value should be scanned and counted. 
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For example, in a typical population of cells, a large 
fraction of those cells are in the resting state (G ) in 
which all of the cells have the same DNA value. In the 
ideal population measurement, a substantial number of 
cells would all have the same value for this particular 
parameter. Using this ideal as a goal, the optical data is 
calibrated by the illumination normalization factor to ob- 
tain a new set of data with a distribution as close to this 
ideal result as possible. A normalization factor can be 
developed and applied to the present sample, to subse- 
quent samples, or, each new sample run can be correct- 
ed based on a normalization factor developed from all 
prior runs. 

In particular, this calibration method functions as fol- 
lows (as illustrated in Fig. 9): As the cells of the sample 
are scanned, a data list is generated for each cell includ- 
ing, jnteralia, the optical neighborhood value (step 300) 
and the location of the maximum value for a select pa- 
rameter (step 302), which represents the center of the 
cell along the scan line, for one or more select parame- 
ters. The neighborhood is then recentered with respect 
to this maximum value (step 304). In this embodiment, 
the optical neighborhood value is not calibrated for other 
factors before being multiplied by the ilium in ization nor- 
malization factor. Next, the optical neighborhood value 
occurring with the highest frequency for the complete 
sample run is determined and designated as the "mode 
optical value" (step 306). Next, the optical neighborhood 
values are reviewed and all cells having optical neigh- 
borhood values within a predetermined range of plus or 
minus "D" of the mode optical value are selected to gen- 
erate a selected subpopulation of cells (step 308). This 
subpopulation of cells is then sorted as a function of lo- 
cation, or pixel position, of the maximum value, or cent- 
er, of each cell along the scan line (step 31 0). The value 
of D is user defined and is typically set to between 2 and 
6 percent of the mode optical value. If this range is se- 
lected to be too broad, it will select cells that do not have 
the constant value, and if this range is too small, then 
not enough cells to create an accurate average will be 
selected. When calibration particles that are known to 
be uniform are used, there is no need to select a sub- 
population of particles within a predetermined range 
around the mode optical calibration value, because ail 
the particles are the same. 

Next, the average optical value of the selected sub- 
population of cells for each pixel position along the scan 
line is determined (step 312). Thereafter, an array of il- 
lumination correction coefficients is computed for each 
scan line including one correction coefficient for each 
location, or pixel position along a scan line (e.g., for 18 
pixels that are preferably used along a scan length L of 
30 samples, there would be 18 coefficients) (step 314). 
Each correction coefficient in the array equals the ratio 
of the corresponding mode optical value to the average 
optical value for the associated parameter. This array is 
the illumination normalization factor. 

The array of correction coefficients is one dimen- 



sional in that it corrects for variations in light excitation 
or collection along the scan line. However, the method 
may also be used to generate such an array for two- 
dimensional systems. For example, the laser beam may 

5 be incrementally moved after each scan over some 
range in a direction perpendicular to the present scan 
beam, rather than moving the microscope stage. At the 
end of this range, the beam movement would be reset 
to zero and the stage incremented. In this case, the re- 

io suiting optical data would be corrected as determined 
by peak location, i.e., the center of the cell in two dimen- 
sions. 

All summed values, i.e., the optical neighborhood 
values, for each parameter are then calibrated for dif- 
15 ferences in illumination along each scan line by finding 
the peak pixel position for each cell and multiplying the 
optical neighborhood value by the normalization factor 
for that scan line position for each parameter (step 316). 
This generates a calibrated set of optical data. 
20 it is desirable to repeat this process for multiple it- 
erations until the data stabilizes. This is accomplished 
by using the calibrated optical data generated by the first 
run as the basis for second and subsequent calibrations. 
In other words, the calibration method is applied to the 
25 initially calibrated optical data stored in the computer 
memory to generate a new and more accurate illumina- 
tion normalization factor, which is then multiplied by the 
initially calibrated optical data to generate a second set 
of more accurately calibrated optical data. This process 
30 can be repeated many times, although the data should 
stabilize after approximately three iterations. The illumi- 
nation normalization factor can be applied directly to 
subsequent sample run optical data or the factor can be 
updated by applying the method to each run and aver- 
ts aging the latest set of correction coefficients from the 
coefficients of prior runs and applying the averaged set 
to the next run. 



40 



Optical Signal Data Acquisition 



After selecting a Protocol and calibrating the instru- 
ment, the operator may initiate one or more data acqui- 
sition runs. The purpose of these runs is to scan an area 
of a test slide defined in the Protocol, find all of the cells 

45 on the slide that meet a threshold criteria for a given 
parameter that is set in the Protocol by the operator, and 
generate a list for each cell found that contains the cor- 
rected optical and morphology values for each param- 
eter of the Protocol, the measurement time, and the X 

50 and Y position of each cell with respect to Home. The 
user places the test specimen on the stage and initiates 
a run by typing a key on the computer keyboard. The 
flow chart of Fig. 3 illustrates the general mechanical 
optical signal data acquisition loop as described above. 

55 The program causes the stage to be driven to Home 
and then moves the stage over the test area. As the 
stage moves, the beam is scanned back and forth to 
create a scan path as shown in Fig. 4. For each X direc- 
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tron scan strip, the optical signal value of each parame- 
ter provided by the specimen is digitized by an analog- 
to-digital (A/D) converter at a sample rate set by the Pro- 
tocol to create a sequence of digital data sample values. 
Typically 120,000 such optical signal values are digi- 
tized and stored in a buffer memory within 3 seconds. 
Typically each scan contains 20 samples so that 2000 
to 6000 forward and reverse scans corresponding to two 
to four parameters with steps of 5 to 10 microns cover 
a length of the same order of magnitude as a slide cov- 
erslip. Data acquisition and processing may be done se- 
quentially or simultaneously. 

The Relationship Between Spot Size and Sampling 
Rate 

In the described embodiment, the optical system is 
designed so that the spot size of the scanning laser 
beam spot 12c at the specimen surface 28 is approxi- 
mately 10 microns when a 20X microscope objective is 
used in the revolving nosepiece. The spot size can be 
changed by rotating the nosepiece to 40X to produce a 
5 micron spot or to 1 0X to produce a 20 micron spot. 

As the beam 12b scans across surface 28, a con- 
tinuous optical analog data stream is produced by each 
of the instrument's optical sensors 20, 22, 24, 26. Figs. 
5a through 5d show graphically how the analog optical 
data frm a fluorescence sensor is sampled by an A/D 
converter. Fig. 5a shows a hypothetical profile of a flu- 
orescence signal from a cell with a level C representing 
a base level of a constituent and with levels representing 
constituent variations C t and C 2 . Fig. 5b shows a scan 
spot intensity profile. The "size - of the spot is measured 
at the half-power points of the Gaussian profile. Fig. 5c 
shows a digital representation of the fluorescence signal 
with inadequate sampling, and Fig. 5d shows a digital 
representation of the fluorescence signal with sampling 
according to the invention. The arrows represent the 
sampling interval. 

The user can define the rate at which each of these 
continuous optical analog data streams are sampled 
and converted by the A/D converter and stored as digital 
data values in the computer memory. The rate can be 
set using the hardware protocol up to a maximum rate 
of 100,000 samples per second. A typical experiment 
uses a rate of 32,000 dual parameter samples per sec- 
ond, i.e., the two analog data streams are each sampled 
at 31 .25 microseconds. 

The velocity of scan beam 12b is fixed by the reso- 
nant frequency of the resonant scanner 70. Such scan- 
ners are commercially available with frequencies to 
2400 Hz. One embodiment uses an 800 Hz. scanner, i. 
e., scan beam 12b moves forward and back, each with 
a distance "L\ the scan length, 800 times per second, 
taking 625 microseconds for the beam to move from one 
end to the other of the scan length L with a cosine ve- 
locity profile. (625nsec/31.25nsec = 20 the number of 
sample points per scan length.) The scan path of the 



laser beam is shown in Fig. 4. 

As shown in Fig. 5b the intensity profile of the scan 
spot 12b on the specimen surface 28 is Gaussian. The 
analog data generated by any sensor from a continuous 
s scan of an optical signal from the specimen is the con- 
volution of the Gaussian intensity profile and the distri- 
bution of the desired cellular constituent signal, i.e., the 
two distributions are multiplied as one is shifted in time. 
Consider a case in which a fluorescence dyed cell 
10 is scanned. A hypothetical profile of fluorescence is 
shown in Fig. 5a which represents levels of fluorescence 
C due to the cell's base level of a constituent and fluo- 
rescence due to constituent variations and C 2 . If the 
sampling rate is low so that the beam moves more than 
75 a spot width between the locations at which the analog 
data stream is sampled, some of the information in the 
fluorescence signal will be lost The digital samples will 
be added later to try to determine the amount of dye 
along a line segment. In the case of low sampling, i.e., 
samples are taken too infrequently, the resultant set of 
values would not be representative of the convolution of 
the distributions of the beam intensity profile and the dye 
distribution because the sampling rate is too low pro- 
ducing "gaps" in the digital data corresponding to the 
analog data as shown in Fig. 5c. Summing these sam- 
ples would underrepresent parts of the base level C and 
furthermore, the information representing the actual in- 
crease due to constituent variation Cz is lost between 
the two peaks representing base level fluorescence C. 
If the sampling rate is such that successive samples in 
the specimen plane are spaced apart by the spot size 
or less, all parts of the specimen will be equally repre- 
sented as shown in Fig. 5d, because Gaussian profile 
spots spaced one spot size apart when added produce 
a uniform distribution along the scan. This is because 
the Gaussian profiles overlap enough when spaced one 
"spot size - apart that, when added, the total equals the 
maximum or peak value of the Gaussian curve. The cell 
specimen is sampled with close enough spot spacing 
so that when the samples from a fluorescence sensor 
are added, with proper calibration, it will yield the accu- 
rate amount of the cells' constituents independent of the 
distribution of that constituent in the cell or the location 
of the cell on the slide. 

The described embodiment uses a conventional cir- 
cuit card available from Data Translation (DT) for con- 
verting analog data to digital values. This card utilizes 
direct memory access (dma) in which the DT card is set 
up by the program to continuously sample each of up to 
four analog signals at a predetermined rate. This DT 
card also converts each of these analog samples to its 
digital data value equivalent, and to store each set of 
signal samples in successive memory locations with a 
predetermined starting memory location, completing the 
operation when a predetermined number of memory lo- 
cations (buffer) is filled with digital values. 

The program sets up this analog-to-digital operation 
by specifying the signals to be sampled, the sampling 
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rate, the initial memory location, and the size of the buff- 
er. One complete X strip is processed from raw data into 
list data in one batch so that the buffer size must be 
equal to the number of scans in a strip (typically 2500) 
times the number of samples in a scan (typically 20) 
times the number of optical signals used in the experi- 
ment (typically 2 to 4); that is, the buffer size is approx- 
imately 100,000. When the buffer is full, program control 
is transferred to the functions that will find cells, find their 
peak values, and determine their optical and morpho- 
logical property values. At this point, there is a buffer 
with typically on the order of 100,000 sequential values 
in it. The dma operation and list value determination op- 
erations can also be overlapped. 

Digital Data Processing 

Once the digital data sample values are stored in 
the computer memory (buffer) the data is processed by 
a variety of protocol controlled functions to correct the 
data, e.g., for background, to calibrate the data, and to 
generate the desired optical and morphological property 
values. These function steps are shown in the flow chart 
of Fig. 7. 

The first data processing program function step, 
200 in Fig. 7, locates the beginning of each forward scan 
(typically 2500 per strip) in.the strip. As shown in Fig. 6, 
one of the signals, typically the forward scatter signal, 
has added to it a synchronization pulse of sufficient am- 
plitude, duration, and negative polarity so that it can be 
distinguished from all normal signals. The synchroniza- 
tion pulse is derived directly from the mirror scan driver 
70, occurring once for each forward scan at the same 
time for every scan and near its beginning so that this 
pulse does not interfere with the actual data. In the de- 
scribed embodiment, only the middle (approximately 
1 2) sample points of the 20 sample scan (i.e., there are 
20 samples per scan length L) are used and will be re- 
ferred to as POINTS. A 30 sample scan length would 
have approximately 18 POINTS. The program first 
searches successive memory locations of the buffer for 
sequences of three values that match the shape of the 
synchronization pulse and produces a table of pointers 
to all of these locations. The starting location of every 
data POINT value is fixed at a known displacement in 
the buffer from each of these pointer locations. Because 
the relationship of the synchronization signal and scan 
position is fixed, the start position of every scan can be 
marked and the data buffer values can each be associ- 
ated with a specific scan position by appropriate record 
keeping in the program. 

Groups of contiguous points for each parameter are 
located in memory based upon their distance from the 
synchronization signals. If parameter samples are taken 
sequentially, i.e., POINT 1 for parameter A, POINT 1 for 
parameter B, POINT 2 for parameter A, etc., then these 
groups of contiguous points will include every other 
sample in memory (for 2 parameters). 



The next function step, 202 in Fig. 7, called by the 
program is responsible for finding cells. The event of lo- 
cating an individual cell is called a "cell find - . The func- 
tion causes a search of each of the POINTS of each 
5 scan to find values for any parameter (with its back- 
ground value subtracted) that exceed the Protocol 
threshold value for that parameter. Thus, any one pa- 
rameter can be used to find a cell by an increase in that 
given parameter's signal value caused by that cell, 
f 0 Background values can be obtained from a calibra- 
tion area as described above or they can be generated, 
either by averaging or determining the minimum values, 
of samples in the same position of scans in one or more 
prior or subsequent neighborhoods. Initially, the pro- 
fs gram computes a table of initial background values for 
the POINTS of a part of a scan strip, i.e., the first ap- 
proximately 50 scans, to determine the lowest value for 
all POINTS for each-parameter signal. The initial back- 
ground value for these POINTS is stored in a table. 

One of the Protocol parameters is the size of the 
neighborhood. The neighborhood has a value N that is 
typically near five POINTS but is adjusted by the user 
to the size and concentration of cells in the test. The 
neighborhood is set to statistically include only one cell. 

In the cell find function, a separate step 204 in Fig. 
7, defines a neighborhood around each cell found. The 
neighborhood value N defines the area (POINTS x 
scans) surrounding the cell find in which other cell finds 
are prohibited, that is, once a single cell is found in a 
given neighborhood, the next find must be outside that 
neighborhood. The number of scans skipped after the 
end of a neighborhood is typically 3. Also, because data 
values are added only within the neighborhoods, cell 
finds are not allowed in areas of the sample at the scan 
edges which comprise less than a full neighborhood. 
The present real time, including date, is taken from the 
computer clock and stored in the cell list. The value of 
N, which is equal to a number of sample POINTS of a 
neighborhood and the number of successive scans of a 
neighborhood is defined in the Protocol. N is preferably 
odd for ease of defining a center of the neighborhood. 

The next function step, 206, compares all values 
along N POINTS of N forward scans, i.e., within one 
neighborhood, centered at the location of the cell find to 
determine the maximum value of that POINT minus its 
corresponding POINT background value for one of the 
Protocol determined signal parameters. The maximum 
value defines the center of the cell and the new location 
of the cell find. The neighborhood is recentered in step 
207. The number of the scan, the POINT number, and 
the strip number are all stored in the main cell list file to 
define this new cell location, along with the computer 
system time. 

In the next step, 208, a new set of POINT back- 
ground values for each neighborhood is computed by 
determining the minimum values of samples for each 
parameter signal using a set number of scans (typically 
5) that end outside the neighborhood, e.g., just prior to 
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the beginning of the current neighborhood. Subtracting 
this background from the digital data POI NTS generates 
optical data POINTS in the computer memory. 

Intensity, velocity and background as determined by 
the calibration run and by the program are all taken into 
account in the next step 210 by multiplying each POINT 
value for each parameter by an array of POINT intensity 
correction factors and velocity correction factors. The 
next function step, 212 in Fig. 7, sums all corrected op- 
tical data values for each parameter for each sample of 
the defined neighborhood centered at the position of the 
peak value for that cell. This sum is the optical neigh- 
borhood value. In the self-calibration method described 
above, this optica I neighborhood value is not initially cal- 
ibrated for intensity or velocity, but is corrected in a sub- 
sequent step. 

In this embodiment, the number of samples of the 
cell taken in each direction is greater than the projection 
of cell length in that dimension divided by the spot size 
in that direction so that when these samples are added, 
the result accurately represents the integrated value of 
that cell parameter and can be directly related to the ab- 
solute quantity of a given cell constituent. This results 
from the principle that the samples are taken at a dis- 
tance apart that is equal to or less than the spot size of 
the incident light beam. The spot size is selected by 
choice of microscope objective power, whereas the 
sampling rate, the step size, and the size of the neigh- 
borhood are all selected in the Protocol to adjust the test 
to the speed, accuracy and sizes of cells tested in the 
specific test. 

The program in the next step, 21 4, generates a dig- 
itized binary neighborhood pattern or morphology word 
for each cell describing its shape. Each POINT in the 
neighborhood for each parameter is compared with that 
parameter's "morphology threshold" set in the Protocol. 
For each POINT where the value exceeds the morphol- 
ogy threshold, a specific position of the morphology 
word in memory is set from zero to one and the word is 
shifted by one. If the threshold is not exceeded, the lo- 
cation is kept at zero and the word is shifted by one. 
Thus, each parameter of each cell is used to generate 
a word equal in length to the area of the neighborhood 
(POINTS x scans) that contains a digital image of the 
cell. In practice, the size of this word can be made small- 
er than the number of POINTS to reduce storage space. 
In a final step, 216, these words are stored in a list for 
each cell along with data for the neighborhood value 
generated in step 21 0, the time of the scan, and cell po- 
sition data. 

Display of Optical and Morphology Parameters 

The above set of functions is repeated for each cell 
find in the buffer representing one X strip. The optical 
property values (generated from the optical neighbor- 
hood value) or the morphological property values (gen- 
erated from the binary neighborhood pattern) of two pa- 



rameters are displayed on a monitor screen as a dot 
whose x position is proportional to one parameter and 
whose y position is proportional to a second parameter. 
When a specific number K of cells are found (typically 
5 500) or a specific number of Y-steps (M) is reached, the 
"run" is completed and a list of values for each cell is 
written in a file to a fixed disk of the computer. The device 
also may be programmed to repeat a run, or rescan, for 
any number P of passes. 
10 jo repeat a run, the microscope stage is moved to 
begin scanning the first X strip and all of the above steps 
are repeated. Simultaneously, other program functions 
control the position of the stage. As in the first pass, the 
stage is initially moved to a HOME position, then to an 
*5 initial start position and ramped up to scanning speed 
while keeping track of position, moved in synchronized 
steps across a strip in the x direction, ramped down in 
speed to a stop, and stepped in the y direction as de- 
scribed above. In this way, by repeat scanning of the 
20 same specimen of cells, the operator can observe kinet- 
ic changes in the individual cells. 

The list of data generated as described above can 
be processed by a data display program that can be run 
at a later time or between strips of the above data ac- 
2S quisition program. The principal features of this program 
are as follows. 

Using a "Data Protocol" the user can define a list of 
properties that are computed by applying mathematical 
functions to one or two optical neighborhood values, 
30 morphology words or time parameters. For example, an 
optical property can be an optical neighborhood value 
for one parameter (e.g., a value that can be converted 
into a number of picograms of DNA), the ratio of two 
such values, the time, or the sum of all of the 1 bits of a 
35 morphology word, i.e., binary neighborhood pattern (the 
cell area). These properties are used as the basis for a 
sequence of cell selections to count and display cell sub- 
populations on the monitor screen. 

Screens are generated on the monitor with axes 
40 corresponding to each of two selected properties. Each 
cell is displayed on the screen as a dot whose position 
on the axis is proportional to a property and whose color 
is proportional to the number of cells occurring at that 
screen position. Using a user controlled mouse, for ex- 
45 ample, a user can bound a selected area (this can be 
extended to additional areas) of the screen by a poly- 
gon. Two additional properties that may be the same or 
different from those of the first screen are designated in 
the data protocol. Cells whose properties are within the 
50 area of the polygon of the first screen are counted and 
displayed on a second screen whose coordinates are 
the second set of properties. A number of polygons, e. 
g., six, can be drawn on the second screen and cells 
within each of these regions are counted. As shown in 
55 Fjg. 8, the counts appear on the screen next to their cor- 
responding regions. 

In addition to counting subpopulations of ceils 
based on successive counts of derived properties, the 
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user can display selected data in a variety of other for- 
mats including isometric projections or contour maps of 
any set of properties. Another aspect of the instrument 
as illustrated in Fig. 8, is the ability of this program to 
generate a polygon on the monitor screen that repre- 
sents a specific location or set of parameter values. The 
program can then instruct the instrument to locate each 
cell whose properties fall within the polygon area. That 
is, each polygon contains a subpopulat ion of cells within 
the sample that exhibits a specific characteristic com- 
mon to each of the cells within the polygon. As each cell 
is automatically positioned under a microscope objec- 
tive 39 for visual observation, a dot on the screen cor- 
responding to that cell within the population is marked 
by a blinking circle centered on it (as shown in the lower 
right hand polygon in Fig. 8). 

After completion of each X strip, the values for each 
cell's parameters are added to a monitor screen display. 
During each test, the user sees each cell represented 
as a dot on the screen at a position determined by the 
values of two parameters or as two such patterns if more 
than two parameters are used. The axes of these graphs 
are labeled using the Protocol parameter names. The 
positions of found cells may also be represented by dots 
on an additional graph on the screen. After a fixed 
number of cells is found, typically about 500, the data 
lists for those cells are stored in a disk file containing a 
header with the Protocol. After either the total area has 
been scanned or the Protocol maximum cell count is 
reached, the remaining cell data is stored in a list file 
followed by the total cell count. 

The user can call for the program to repeat the scan 
one or more times. In these subsequent passes, all of 
the above functions, with the optional exception of cell 
finding, are repeated. The locations of all found cells 
from the first pass are read into a table in memory and 
these locations are used at the appropriate times to in- 
dicate the starting position of the peak value finding rou- 
tine. The new time is stored in the new list and all other 
functions are identical to the first pass. 

Data Analysis 

The result of the program is a file containing a set 
of data for typically 1000 to 100,000 cells. These data 
values are used as the input to an analysis and display 
program which can be identical to any of the programs 
described in the flow cytometry literature. These pro- 
grams are reviewed in the texts, Melamed, et al. Flow 
Cytometry and Sorting (John Wiley & Sons, 1979, Ch. 
20), or Shapiro, H., Practical Flow Cytometry, (A.K. Liss, 
1985, Ch. 5). " 

The monitor screen displays describing morpholo- 
gy, time and position, however, are another feature. For 
example, time can be represented graphically as one 
additional parameter along with scatter or fluorescence 
parameters or specific functions of parameters may be 
tested against time variable thresholds or displayed as 



a function of time to implement experiments or tests. 
Time can also be displayed by modifying dot or contour 
colors on a display to mark the passage of time. Func- 
tions of the morphology word such as the total one bit 

5 count indicative of area of that parameter may be used 
like any parameter. It is contemplated that the large 
number of total parameters that can be obtained from 
the morphology and value data can be best analyzed 
with techniques described in the pattern recognition or 

io artificial intelligence literature to achieve the end goal of 
identifying specific cell types or pathological cells on a 
slide. 

Once cells with unique parameter values are found 
by the analysis program, the test specimen can be 
placed on the microscope stage and a program called 
to automatically center the stage on each of these 
unique cells in turn by reading in a list of the X and Y 
values of cells to be reviewed. Properties of these cells 
may be displayed on the monitor screen. The user may 
use this to verify the automatic data analysis, examine 
the cells to further determine their pathology, check the 
quality of the manual procedures in a clinical laboratory, 
or physically separate or sort cells of interest using cook- 
ie cutter like devices. 

The kinetic behavior of live cells, primarily blood 
lymphocytes, when reacted with certain antigens during 
a test, may be an important new method of testing for 
systemic pathologies such as autoimmune disease or 
cancer. As described above, cells with specific proper- 
ties may be identified and their locations stored in mem- 
ory. The specimen in a cuvette or chamber on the stage 
can be rescanned and the new values of these selected 
cells recalculated and listed. 

The device also contemplates sequential analysis 
techniques to maximize speed of analysis. Cells of in- 
terest found during a first pass will be examined for other 
parameters or at greater resolutions to determine prop- 
erties that require greater analysis complexity and time. 
Other embodiments are feasible. 



Claims 

1. A method for characterizing a population of cells, 
the method comprising: 

(a) scanning the cell population on a surface 
along a scan line with a laser beam with a spot 
size comparable in size to the cells to be 
scanned to produce analog signals from a sen- 
sor that detects light emitted by the cells and 
sampling said analog signals to produce digital 
data samples, said samples being taken at a 
sampling rate set such that the distance trav- 
elled per time as the spot is scanned between 
consecutive sampled locations approximates 
said spot size; 

(b) displacing the scan line of said laser with 
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respect to said surface such that different said 
samples represent different locations within the 
cell population; and 

(c) processing the digital data to generate cell 
constituent data that characterizes the cells in s 
the population. 

A method according to Claim 1 , wherein said spot 
size is 3 to 50 microns. 

10 

A method according to Claim 1, wherein said 
processing step comprises correcting said digital 
data samples with a background level to generate 
said optical data. 

15 

An apparatus for measuring and displaying cellular 
properties of a population of cells comprising: 

a laser source for generating a beam of light to 
illuminate the cells to be measured with a beam 20 
spot comparable in size to the cells to be 
scanned, an optical signal being emitted as a 
result of said spot illuminating said cells; 
a surface upon which the cells are located; 
an optical path for directing said spot from said 25 
beam source to the cells on said surface; 
a scanner interposed in said optical path be- 
tween said beam source and the cells for scan- 
ning said spot across said surface along a scan 
line; 30 
a sensor for measuring said optical signal; 
an analog-to-digital converter arranged to sam- 
ple said optical signal at a specified rate and to 
produce digital data; whereby said sampling 
rate is set such that the distance travelled per 35 
time as the spot is scanned between consecu- 
tive sampled locations approximates said spot 
size; 

means for displacing the scan line of said 40 
laser with respect to said surface such that 
said laser spot scans continuously over a 
plurality of cells on said surface; 
a data processor for deriving cell constitu- 
ent data from said digital data correspond- 45 
ing to individual cells, and processing said 
cell constituent data to display cellular 
properties of the ceil population; and 
a memory for storing said digital data pro- 
duced by said converter. so 

An apparatus according to Claim 4, wherein said 
spot size is 3 to 50 microns. 

An apparatus according to Claim 4, wherein said 55 
scanner is a resonant galvanometer scanner. 

An apparatus according to Claim 4, wherein the 



wavelength distribution of said beam is controlled 
by the data processor. 



Patentanspruche 

1 . Verfahren zur Charakterisierung einer Zellpopulati- 
on, das folgendes umfaBt: 

(a) Scannen der Ze I {population mit einem La- 
serstrahl mit einer FleckgrdGe, die mit der Gro- 
Be der zu scannenden Zellen ubereinstimmt, 
auf einer Oberflache entlang einer Scanlinie 
zur Erzeugung von analogen Signalen durch 
einen Sensor, der zur Messung des von den 
Zellen abgestrahtten Lichts dient, und Abtasten 
der analgogen Signals fur die Erzeugung von 
Abtastwerten in Form von digitalen Daten, wo- 
bei die Abtastwerte mit einer Abtastrate gewon- 
nen werden, die so eingesteilt ist, daft der in 
dem Zeitraum zum Abtasten zweier aufeinan- 
derfolgender abgetasteter Stellen durch den 
Fleck zuruckgelegte Abstand in etwa der Fleck- 
groBe entspricht; 

(b) Verschieben der Scanlinie des Lasers auf 
der Oberflache in einer solchen Weise, daB 
verschiedene Abtastwerte verschiedene Stel- 
len innerhalbderZellpopulation reprasentieren 
und 

(c) Verarbeiten der digitalen Daten zur Erzeu- 
gung von Daten zu den Zellbestandteilen, die 
die Zellen in der Population kennzeichnen. 

2. Verfahren nach Anspruch 1 , bei dem die Fleckgrd- 
Be 3 bis 50 Mikrometer betragt. 

3. Verfahren nach Anspruch 1 , bei dem zu diesem Ver- 
arbeitungsschritt die Korrektur der in Form von di- 
gitalen Daten vorliegenden Abtastwerte mittels ei- 
nes Hintergrundpegels zur Gewinnung der opti- 
schen Daten gehort. 

4. Apparatur zum Messen und Anzeigen der Zellei- 
genschaften einer Zellpopulation, bestehend aus: 

einer Laserquelle fur die Erzeugung eines 
Lichtstrahls zur Beleuchtung der Zellen, die mit 
Hilfe eines Strahlflecks zu untersuchen sind, 
der in der GroBe vergleichbar ist mit den zu 
scannenden Zellen, wobei infolge der Beleuch- 
tung der Zellen mit Hilfe des Strahlflecks ein op- 
tisches Signal abgegeben wird; 

einer Oberflache, auf der die Zellen angeordnet 
sind; 
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einer Einrichtung zum Verschieben der Scan I i- 
nie des Lasers auf der Oberflache, so daB der 
Laserfleck kontinuierlich auf der Oberflache 
uber eine Anzahl von Zellen gefOhrt wird; 

einer Datenverarbeitungseinheit zum Abieiten 
der zu den Zellbestandteilen gewonnenen Da- 
ten aus den den einzelnen Zellen entsprechen- 
den digitalen Daten und zum Verarbeiten der 
zu den Zellbestandteilen gewonnenen Daten, 
so daft die Zelleigenschaften der Zellpopulati- 
on angezeigt werden konnen und 

einem Speicher zum Speichern der mit Hilfe 
des Wandlers gewonnenen digitalen Daten. 

5. Apparatur nach Anspruch 4, bei der die FleckgroGe 
3 bis 50 Mikrometer betragt. 

6. Apparatur nach Anspruch 4, bei der der Scanner 
ein Resonanz-Galvanometerscanner ist. 

7. Apparatur nach Anspruch 4, bei der die Wellenlan- 
genverteilung des Strahls durch die Datenverarbei- 
tungseinheit gesteuert wird. 



Revendications 

1. Precede pour caract6riser une population de cellu- 
les, le procede comprenant les Stapes suivantes : 

(a) balayage de ia population de cellules sur 
une surface le long d'une ligne de balayage par 
un faisceau laser ayant une dimension du spot 
comparable a la dimension des cellules a ba- 
layer pour former des signaux analogiques au 



moyen d'un capteur qui detecte la lumiere 6mi- 
se par les cellules, et 6chantillonnage desdits 
signaux analogiques afin de former des echan- 
tillons de donnees, lesdits echantillons etant 
preleves a une Vitesse d'echantillonnage telle 
que la distance parcourue pendant le temps de 
balayage du spot entre des emplacements 
d'echantillonnage successes approche ladite 
dimension du spot ; 

(b) defacement de la ligne de balayage dudit 
laser par rapport a ladite surface de maniere 
que lesdits difterents echantillons represented 
differents emplacements dans la population 
des cellules ; et 

(c) traitement des donn6es numeriques pour 
g6ne>er des donn6es de constituants des cel- 
lules qui caracterisent les cellules de la popu- 
lation. 

2. ProcedS selon la revendication 1 , dans lequel ladite 
dimension du spot est comprise entre 3 et 50 mi- 
crometres. 

Proced6 selon la revendication 1 , dans lequel ladite 
6tape de traitement comprend une correction des- 
dits echantillons de donnees numeriques d'un ni- 
veau de fond pour gSnerer lesdites donnees opti- 
ques. 

30 4. Appareil pour mesurer et afficher des propriety cel- 
lulaires d'une population de cellules comprenant : 

une source laser pour gendrer un faisceau lu- 
mineux qui illumine les cellules a mesurer par 
55 un faisceau ayant un spot de dimension com- 

parable a celle des cellules a balayer, un signal 
optique 6tant emis du fait de rillumination des- 
dites cellules par (edit spot ; 
une surface sur laquelle les cellules sont 
40 placees ; 

un chemin optique pour diriger ledit spot emis 
par ladite source de faisceau vers les cellules 
placees sur ladite surface ; 
un scanner intercale sur ledit chemin optique 
45 entre ladite source de faisceau et les cellules 

pour balayer ledit spot a t ravers ladite surface 
le long d'une ligne de balayage ; 
un capteur pour mesurer ledit signal optique ; 
un convertisseur analogique/numerique agen- 
50 ce pour echantillonner ledit signal optique a 

une Vitesse specified et produire des donn6es 
numeriques ; grace a quoi ladite vitesse 
d'echantillonnage est fix6e de facon que la dis- 
tance parcourue par unite de temps quand le 
55 spot est balaye entre deux emplacements 

echantillonn6s successifs approche ladite di- 
mension du spot ; 

un moyen pour deptacer la ligne de balayage 



einem Strahlenveriauf, auf dem der Fleck von 
der Strahlquelle auf die Zellen auf der Oberfla- 
che gelenkt wird; 

einem Scanner im Strahlenveriauf zwischen s 
der Strahlquelle und den Zellen zum FOhren 
des Flecks auf einer Scanlinie uber die Ober- 
flache; 

einem Sensor zum Messen des optischen Si- 10 
gnals; 

einem Analog-Digital-Wandler zum Abtasten 
des optischen Signals mit einer bestimmten 
Abtastrate undzur Erzeugung von digitalen Da- '5 
ten, wobei die Abtastrate so eingestellt wird, 
daB der in dem Zeitraum zum Abtasten zweier 
aufeinanderfolgender abgetasteter Stellen 
durch den Fleck zuruckgelegte Abstand in etwa 
der FleckgrdGe entspricht; 20 
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dudit laser par rapport a ladite surface, de ma- 
niere que ledit spot laser effectue un balayage 
continu de la plurality des cellules placees sur 
ladite surface ; 

un processeur de donn6es pour dSduire des 5 
donnees de constituants des cellules desdites 
donnees num6riques correspondant aux cellu- 
les indMduelles, et pour traiter lesdites don- 
nees de constituants des cellules afin d'afficher 
des proprietes cellulaires de la population de 10 
cellules ; et 

une m6moire pour memoriser lesdites donnees 
numeriques produites par ledit convertisseur. 

Appareil selon la revendication 4, dans lequel ladite is 
dimension du spot est comprise entre 3 et 50 mi- 
crometres. 

Appareil selon la revendication 4, dans lequel ledit 
scanner est un scanner a galvanometre resonnant. 20 

Appareil selon la revendication 4, dans lequel la dis- 
tribution des longueurs d'ondes dudit faisceau est 
commandee par le processeur de donnees. 
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Abstract of EP0421736 

The invention relates to an apparatus and 
methods for generating multiparameter optical 
data that characterize a population of cells. The 
invention includes the steps of scanning the cell 
population with a beam to produce sets of digital 
data samples, each sample set of digital data 
representing multiparameter optical interactions 
from a specific location within the cell population; 
storing the digital data, e.g., in a computer 
memory; locating a cell within the population, 
e.g., by comparing the digital data to a 
preselected threshold value; defining a 
neighborhood around the digital data 
representing the located cell; estimating a 
background level for the neighborhood based 
upon digital data corresponding to locations 
outside the neighborhood; and correcting each of 
the samples corresponding to the neighborhood 
with the estimated neighborhood background 
level to generate the optical data. The invention 
further relates to specific methods of background 
correction and data calibration as well as specific 
sampling features to enable precise estimates of 
multiple cellular constituents or other cell 
properties at high rates of speed. 
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